Quantum fidelity for one-dimensional Dirac fermions and two-dimensional Kitaev 

model in the thermodynamic limit 
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We study the scaling behavior of the fidehty (F) in the thermodynamic hmit using the examples 
of a system of Dirac fermions in one dimension and the Kitaev model on a honeycomb lattice. We 
show that the thermodynamic fidelity inside the gapless as well as gapped phases follow power- 
law scalings, with the power given by some of the critical exponents of the system. The generic 
scaling forms of F for an anisotropic quantum critical point for both thermodynamic and non- 
thermodynamic limits have been derived and verified for the Kitaev model. The interesting scaling 
behavior of F inside the gapless phase of the Kitaev model is also discussed. Finally, we consider 
a rotation of each spin in the Kitaev model around the z axis and calculate F through the overlap 
between the ground states for angle of rotation 77 and 77 + dr], respectively. We thereby show that 
the associated geometric phase vanishes. We have supplemented our analytical calculations with 
numerical simulations wherever necessary. 
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I. INTRODUCTION 

A quantum phase transitionir— driven exclusively by 
quantum fluctuations at zero temperature is associated 
with a dramatic change in the symmetry of the ground 
state of a many-body quantum Hamiltonian. A num- 
ber of measures from quantum information theory such 
as entanglement^, entanglement entropy^i^, Loschmidt 
echoiS, decoherenceii and quantum discordi^ii^ are able 
to capture the singularities associated with a quantum 
critical point (QCP). Consequently, there is a recent 
upsurge in the investigation of quantum critical sys- 
tems from the perspective of quantum information the- 
ory in an attempt to establish a bridge between these two 
fieldsii-i^. 

An important information theoretic concept that is be- 
ing investigated extensively for quantum critical systems 
is the quantum fidelity (F)i2r— . Let us consider a d- 
dimcnsional Hamiltonian H{X) which contains an exter- 
nally tunable parameter A, such that the system is at a 
QCP when A = 0. Considering two ground state wave 
functions \ipoW) ^^nd |'0o(-^ + <5))i which are infinitesi- 
mally separated in the parameter space as 6 — > 0, wc 
define the fidelity as 



FiX,X + S) = \{M>^)\M^ + 5))\ 



(1) 



where L is the linear dimension of the system, and S de- 
notes a small change in the parameter A. The first non- 
vanishing term in the expansion of F, namely, the fidelity 
susceptibility provides a quantitative measure of the 
rate of change of the ground state under an infinitesi- 
mal variation of A. By exhibiting a sharp decay around 
a QCP even for a finite size system, the fidelity turns 
out to be one of the fundamental probes for detecting 
the ground state singularities associated with a quantum 



phase transition without making reference to an order pa- 
rameter. At the same time, the fidelity susceptibility de- 
fined through the relation xf = -i2/L'^){hiF)/6^\s^a = 
— {\/L'^){d'^F/d5'^) usually diverges with the system size 
in a universal power-law fashion with an exponent given 
in terms of some of the quantum critical exponents. For 
a marginal or relevant perturbation A, the application of 
the adiabatic perturbation theory leads to a generic scal- 
ing formi^>26i27^^ Qf gjygj^ Y>y xf{X = 0) - L'^/^-'^ 
at the QCP, whereas away from the QCP, the scaling 
changes to xf ~ |A|'^'*~^ for L > A^"; here is the 
critical exponent describing the divergence of the spatial 
correlation length close to the QCP, i.e., ^ ~ . 

While the fidelity susceptibility approach usually as- 
sumes small L and (5 — >■ 0, the fidelity per site^^"— has 
been calculated in the thermodynamic limit for an ar- 
bitrary value of S in some recent studies^^"— . In such 
cases, the fidelity differs significantly from unity unlike 
in Eq. ([T]). The fidelity per site is also able to indicate 
the appearance of a quantum phase transition in the spin- 
1/2 XY model in a transverse magnetic field. Another 
measure of fidelity applied to mixed states, namely, the 
reduced fidelity^"~ has also provided important insights 
into quantum critical phenomena. We note that similar 
studies have been carried out on the scaling of the geo- 
metric phase^i^ which is closely related to the fidelity 
susceptibilitjii^ close to critical^"— and multicritical^ 
points. 

As seen from Eq. ^ and also the preceding discus- 
sion, the knowledge of xf should be sufficient to draw 
conclusions about the behavior of F and the associated 
QCP for small system sizes and in the limit S ^ 
[S'^L'^Xf/'^ ^ !)■ However, in the thermodynamic limit 
(L — > 00 at fixed S), the expansion in Eq. ([l} up to the 
lowest order becomes insufficient. Recently, Rams and 
Damski^ have proposed a generic scaling relation^^ valid 
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in the thermodynamic hmit given by 



II. DIRAC FERMIONS IN ONE DIMENSION 



\nF{X- 5,X + d) ~ -L'^\5\'"^A 



(2) 



where A is a scahng function; this relation interpolates 
between the fidelity susceptibility approach and the fi- 
delity per site approach. In deriving the scaling relation 
in Eq. (I2|), it is assumed that the fidelity per site is well 
behaved in the limit L — > oo. the QCP is determined 
by a single set of critical exponents, and i/d > 2 so that 
non-universal corrections are subleading^^. In particu- 
lar, at the critical point A = 0, the fidelity, measured 
between the ground states at +S and —5, is non-analytic 
in S and satisfies the scaling InF —L'^\S\'^'^. On the 
other hand, away from the QCP, i.e., for \S\ <^ |A| ^ I, 
the scaling gets modified to Ini^ ^ ~L'^P\X\'^'^~^ . This 
scaling has been verified for an isolated quantum crit- 
ical point using one-dimensional transverse Ising and 
XY Hamiltonians^. Moreover, near a QCP a cross- 
over has been observed from the thermodynamic limit 
(i|(5|'^ ^ 1) to the non-thcrmodynamic (small system) 
limit ^ 1) where the concept of fidelity suscepti- 

bility becomes useful. We note that Eq. ([2]) is an example 
of the Anderson orthogonality catastrophe^ which states 
that the overlap of two states vanishes in the thermody- 
namic limit irrespective of their proximity to a QCP. 

In this paper we investigate the scaling of the thermo- 
dynamic fidelity in a one-dimensional system of Dirac 
fcrmions with a mass perturbation^"— and the two- 
dimensional Kitaev model on a honeycomb lattice^i^ 
close to or inside the gapless phases of their phase 
diagrams, thereby extending previous studies to more 
generic situations. For the onc-dimcnsional system, we 
have verified the scaling predicted in Ref. |33|. We also 
propose a generic scaling form for the thermodynamic fi- 
delity in the vicinity of an anisotropic quantum critical 
point (AQCP), and we verify it for the AQCP present in 
the Kitaev model phase diagram^. 

The paper is organized in the following way. In Sec. 
II, we analytically derive the scaling relations of the ther- 
modynamic fidelity for non-interacting spinless massive 
Dirac fermions in one dimension and propose a gener- 
alization to the case of interacting fermions (called a 
Tomonaga-Luttinger liquid). In Sec. Ill, wc concentrate 
on the fidelity of the Kitaev model on the hexagonal lat- 
tice for different values of the coupling parameters. Wc 
derive the scaling laws for both the thermodynamic and 
non-thermodynamic limits for an AQCP and verify these 
numerically for the Kitaev model. In Sec. IV, we calcu- 
late the overlap between two ground states of the Kitaev 
model; in one state, all the spins are rotated about the z 
axis by an angle 77, while in the other, they are rotated by 
an angle rj -\- drj. Wc thus derive the form of the fidelity 
through an expansion in powers of drj. 



In this section we consider a system of spinless Dirac 
fermions in one dimension with a mass perturbation and 
verify the scaling of the thermodynamic fidelity as pre- 
dicted in Ref. |3J. Let us first consider non-interacting 
fermions. The Hamiltonian we consider is 



H = 



k>0 



(cic-fc + clfcCfe) , (3) 



where cj^ (ck) is the fermionic creation (annihilation) op- 
erator for wave vector fc, m is the mass, and we have set 
the velocity v — 1 for convenience. In the two-level sys- 



tem given by c\ck 
the form 

H 



where hk 



; = 1, the Hamiltonian takes 



E ( 4 cU 

k>0 

k m 
m —k 



Ck 
C-k 



(4) 



The normalized ground state of this is given by 

ip{k, m) 



y^2 (fc2 + to2) + 2fcVPT^ 



-yk^ + m? — k J ' 
(5) 



with the energy Ek = —y'k^ + m?. 

If we now consider two systems with masses mi and 
m2, the fidelity between the two ground states is given 
by 

t (nil, 1712 ) = Yl mk,miMik,m2))\. (6) 

fc>0 

Note that we have taken k to be strictly positive in all 
the equations above. It is important to exclude the mode 
with fc = 0, otherwise the fidelity is exactly equal to 
zero if mi and 7712 have opposite signs; this is because 
(-0(0, mi)|7/'(0, 7712)) = if 771177^2 < 0. The simplest way 
to exclude a zero momentum mode is to impose antiperi- 
odic boundary conditions, 'ijj{x = L) = —'ip{x = 0), so 
that kn = {n/ L){2n + 1), where L is the system size and 
77 = 0, 1, 2, • • • . (Note that the spacing between succes- 
sive values of k is 2t:/L). We can then write Eq. © 
as 

00 

i^(77ll, 7772) = |('0(fc„, 777l)|7/;(fc„,77l2))|. (7) 

Now we consider the case with 7?7i = 77i and 7772 = 
—777 so that the states lie on the two sides of the gapless 
critical point, and 777 plays the role of 5 discussed in the 
previous section. We find that 



{'ij}{kn,m)\i}j{kn, -777)) ^ 



kl + kny/kl + 7772 



1,2 



(8) 
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Since fc„ = (7r/i)(2n + 1), we see that the fidelity is a 
function of a single parameter given by mL. We now 
consider two cases: (i) mL ^ 1 and (ii) mL ^1. In 
both cases, we will assume that L ^ 1. (Cases (i) and 
(ii) will be respectively called the thermodynamic and 
non-thermodynamic limits in the next section). 

In case (i), we can take fc to be a continuous variable 
so that the fidelity is given by an integral, 



F(m, 



cxp 



L 



dk 
2^ 



In |(V'(fc, fni)\'ip{k^ 771,2) 



(9) 

By writing k ~ mx in the integral, wc find that 
F{m, —to) is given by e~""^, where 



dx , 
In 

2tt 



+ x\/x^+l 



+ 1 + xVx'^ + 1 



(10) 



We conclude that the fidelity satisfies the scaling form 
\nF ^ —cLm which is in agreement with the predictions^ 
that InF ~cL6'^'^ where v = d = 1 and m = 5 \n the 
present case. 

In case (ii), we can expand (?/'(fc„, TO)|'i/'(fc„, — to)) ~ 
1 — TO^L^/[2(2n + 1)^] to lowest order in mL. Wc then 
obtain 



ln_F'(TO, —to) 



m^L^ 



00 

E 

n=0 



TT^TO^L^ 



(2n + l)2 16 ■ 

(11) 

Hence wc find the scaling relation Ini^ ~ S^L^^'^ in the 
non-thcrmodynamic limit; we can further conclude that 

XF ^ L ' 

We now consider what happens if the fermions were in- 
teracting; in one dimension, such a system is described by 
Tomonaga-Luttinger liquid theory^"—. To be specific, 
let us consider the spin-1/2 XX Z chain in a transverse 
magnetic field; the Hamiltonian is given by 



H 



00 

E 

I— —00 



Ota. 



T Z Z 

'Jz"n"n+1 



(12) 

We first set /i„ = 0. Upon using the Jordan- Wigner 
transformation which takes us from spin-1/2 to spinless 
fermions in one dimension^, we find that the first two 
terms in Eq. (l/2)(cr>^+i + o-^cr^+i), lead to a 

tight-binding Hamiltonian of the form — 4- 1 + 

c|i+iC„); Fourier transforming to fc-space, and linearizing 
around the two Fermi points lying at A: = ±7r/2 gives 
the first term in Eq. (jS]). The third term in Eq. ([T2|) . 
cr^cr^+i, leads to a four-fermion interaction of the form 

c-li'^nCn+iCn+i- Eq- dl) describes a system of massless 
interacting fermions if Jz lies in the range — 1 < Jz < 1. 
(The cases Jz = —1 and 1 describe an isotropic ferro- 
magnet and isotropic antifcrromagnet respectively) . The 
system is characterized by a Luttinger parameter K given 

by 



so that K goes from cxd to 1/2 as goes from —1 to 1. 
If = 0; there are no interactions between the fermions 
and we obtain K = 1. 

We now introduce an alternating magnetic field of the 
form hn = to(— 1)"; this introduces a coupling between 
the modes at the Fermi points k — ±7r/2 and there- 
fore gives rise to the mass term in Eq. Q if = 
{K = 1). The most efficient way of studying the low- 
energy, long wavelength modes (i.e., the modes near 
the Fermi points) of a one-dimensional system of in- 
teracting fermions {K ^ 1) is to use the technique of 
bosonization^"— . Bosonization uses a quantum field 
theory which is defined in terms of a scalar field 0. In 
this description, the action is given by 



2K 



dtdx 



dt 



dx 



(14) 



for TO = (we have again set the velocity equal to 1); 
the effect of interactions appears through the Luttinger 
parameter K. The contribution of the mass term to the 
action takes the form 



dtdx TO cos(2v'7r(/)). 



(15) 



The operator cos(2-\/7r0) is known to have mass dimen- 
sion K] let us denote the coefficient of this operator 
in the action by A. It turns out that A effectively be- 
comes dependent on the length scale L, and A(L) satis- 
fies the renormalization group (RG) equation dA/dln L = 
(2 — K)\. Given the initial value of A(a) = to at some 
microscopic length scale a (such as the lattice spacing), 
and assuming that to ^ 1, the RG equation implies that 
A(i) grows and becomes of order 1 at a length scale given 
by ^, where £,/a ^ l/m^^^'^"-^K Hence the mass gap of 
the theory is given by 1/^ ~ jti^/^^^^^ leading to a low- 
energy dispersion given by Uk — \/fc^~{-TO^7(2^-^. One 
can then argue qualitatively that our above arguments 
about fidelity would remain valid, but with a renormal- 
ized mass term given by m^^^'^~^\ For case (i) where 
Lm^/^'^~^^ » 1, we would eventually find that the fi- 
delity scales as lnF(TO, —to) ~ — c Lm^/'^'^~^\ where c' is 
a prefactor which differs from c due to the presence of the 
interactions. Noting that the correlation length exponent 
in the presence of the mass perturbation \sv= 1/(2 — A'), 
we find that the scaling in Eq. ^ should also hold good, 
with d = \ and S ~ m. The above analysis is only valid 
for K < 2; for K > 2, the mass perturbation is irrelevant 
in the sense of the renormalization group, and the fidelity 
has to be calculated in some other way which we will not 
pursue here. 

We will not examine here the effects of a perturbation 
which takes the system across a Kosterlitz-Thouless tran- 
sition; this occurs if /i„ = and Jz crosses 1. The fidelity 
across this transition is considerably harder to analyze, 
and we refer to the work done by different groups^^— . 
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III. KITAEV MODEL 



In this section, we will exploit the solvability of the Ki- 
taev model on the honeycomb lattice^ to calculate the 
fidelity of the model as a function of various parameters 
in the thermodynamic limit. We note that the fidelity per 
site studied for the same model has been able to detect 
quantum phase transitions^^, and the fidelity susceptibil- 
ity has been calculated previously in the limit of small 
system sizes^iiS^. 



> GAPPED 




> GAPLESS 



A. Model, phase diagram and fidelity 



The Hamiltonian of the Kitaev model is given by 



H 



E 

j-\-l — even 



3,1^3 + 1,1 



(16) 

where j and I respectively denote the column and row 
indices of a honeycomb lattice (see Fig. i)^ ~52,59 ^ 
will assume that the couplings J^, for i = 1,2,3, are all 
positive; if some of them are negative, they can be made 
positive by appropriate tt rotations about the x, y or z 
spin axis. For the moment, our study will be restricted 
to the case Ji = J2, although we will comment on the 
case with Ji ^ J2 in Sec. Ill E. 




7-1 j j+l j+2 ,/+3 j+4 j+5 



FIG. 1: (Color online) Schematic representation of the Ki- 
taev model on a honeycomb lattice showing the bonds with 
couplings Ji, J2 and J3. Mi and M2 are spanning vectors of 
the lattice. Sites 'a^' and '63' represent the two inequivalent 
sites which make up a unit cell. (After reference [s^) 



FIG. 2: (Color online) Phase diagram of the Kitaev model; all 
the points in the triangle satisfy Ji + J2 -f J3 = 1. The gapless 
phase is the region in which the couplings satisfy the triangle 
inequalities given by Ji < J2 + J3, J2 < J3 + Ji and J3 < 
Ji + J2, i.e., the points inside the inner equilateral triangle. 
Along the dashed vertical line J3 is varied holding Ji = J2, 
and the anisotropic quantum critical point (A) at J3 = J3,c = 
Ji + J2 ~ 2Ji is indicated. Our focus is to calculate the 
fidelity between ground states lying on this vertical line. 



We define the Jordan- Wigner transformation as 



aj,i = 



b3,i = 



6' 



\l — — OQ / 

n <i 1 '^h fo'^ *5ven j + l, 

- — 00 / 

li <) <i for odd J + I, 

— — 00 / 

ri for odd J + 1 

——00 / 



\i— — oo 

f J-1 



(17) 



where ttj^;, a^ j, hj^i and 6^-; are all Majorana fermions, 
i.e., they are Hermitian, their square is equal to 1, and 
they anticommute with each other. Instead of using the 
indices (j, X) to specify the sites, we can use the two- 
dimensional vectors n ~ VSi'T-i + (-^^ + \3)'ni which 
denote the midpoints of the vertical bonds of the honey- 
comb lattice; here i denotes the unit vector along the hor- 
izontal (labeled by j, j -l-1, etc) and similarly j is the unit 
vector along the vertical direction. Here n\ and ^2 run 
over all integers so that the vectors n form a triangular 
lattice. The Majorana fermions a^i (a^) and 6^ (6^) are 
located at the bottom and top lattice sites respectively of 
the bond labeled by n. The vectors M\ = -^i — ij and 



A/2 



+■ |j shown in Fig. 1 are the spanning vectors 
of the lattice. (We have set the nearest neighbor lattice 
spacing to unity). 

The Fourier transforms of the Majorana fermions are 



5 



given by 



E 



ik-n 



t —ik-n 

al e 

k 



(18) 



satisfying {aj:,at^} ~ ^^k'^ ^'^'^ similarly for a^, and 
6^. In Eq. (|18p . L is the number of sites (hence the num- 
ber of unit cells is L/2), and the sum over k extends over 
half the Brillouin zone of the hexagonal lattice because of 
the Majorana nature of the fermion o^^i^^ . The full Bril- 
louin zone is given by a rhombus with vertices lying at 
{h^^ky) = (±27^/^/3,0) and (0,±27r/3); half the Brillouin 
zone is given by an equilateral triangle with vertices at 
{k^.ky) = (27r/v^,0) and (0,±27r/3). 

In terms of the Majorana fermions, the Hamiltonian in 
Eq. ([T5|) takes the form 



(19) 

where Dfi = i 6^a^. We note that the operators Dfi 
have eigenvalues ±1, and commute with each other and 
with H'; hence all the eigenstates of H' can be labeled 
by specific values of 0,-%. (We observe that the Hamilto- 
nian H' gives dynamics to the fermions an and 6h, but 
the fermions a'- and 6'- have no dynamics since ib'^a'^ is 
fixed). The ground state can be shown to correspond to 
Df{ = 1 for all fiiM For Df{ = 1, the Hamiltonian can be 
diagonalized into the form 



(20) 



where can be written in terms of Pauli matrices as 

~ ^k '^^ l^k '^'^ ^ 
where = 2[Ji sin(fc • il/i) — J2 sin(A: • Af2)], 
and I3j: = 2[J3 + Ji cos(fc • Mi) + J2 cos(fc • Nh)]- 

(21) 

The energy spectrum of H' consists of two bands with 
energies given by 



± 



(22) 



The energy gap Et — E^ vanishes for specific values of 

k when | Ji — J2 1 < < Ji -\- J2 giving rise to a gapless 
phase of the model. The gapless and gapped phases of 
the model arc shown in Fig. [2] in terms of points in an 
equilateral triangle which satisfy Ji + J2 + = 1 and 
all Ji > 0, with the value of Ji being given by the dis- 
tance from the opposite side of the triangle as indicated 
by the arrows. In the limit J3 = 0, the Hamiltonian 
([T6| reduces to a one-dimensional version of the Kitaev 
model^ which in turn can be mapped to the transverse 
Ising chain following a duality transformatioi*^. 



In particular, let us consider the critical line = Ji + 
J2 which separates one of the gapped phases from the 
gapless phase in Fig. [2j On this line, the energy vanishes 
at the three corners of half the Brillouin zone given by 
k = (27r/-\/3, 0) and (0, ±27r/3); these three points are 
actually equivalent to each other because they are related 
by shifts by the lattice vectors G± = (27r/A/3, ±27r/3). If 
{dkx, dky) denotes a small deviation from any one of these 
three points, we find that the energy is highly anisotropic 
with respect to this deviation. Namely, for J3 = Ji + J2, 
the quantities and /3g appearing in Eqs. (|2m22p are 
given by 



V3(J2 - Ji)dh 



J2)dky 




h = 



respectively, to lowest order in dk^ and dky. We see 
that varies linearly in one particular direction in the 
plane of (dk^^dky), while f3j: varies quadratically in any 
direction. We thus have an AQCB^. For simplicity, 
we will mainly restrict our attention below to the case 
where Ji = J2 is held fixed and J3 is varied along the 
dashed vertical line shown in Fig. [2l Then the point J3 = 
J3_c = 2 Ji marked by A is an AQCP, with the energy gap 
vanishing near the three points as Ej: {dk^)^ and dky 
for deviations along the k^ and ky directions respectively. 
For Ji — the dispersion is linear along the vertical 
direction j and quadratic along the horizontal direction 
i (see Fig. [T]). This implies, for the analysis given in 
Sec. HI B, that the correlation length exponent v± = 1 
and Lj^ is the length of the system in the j direction, 
while the exponent v\\ = 1/2 and L\\ is the length of 
the system in the i direction. For a more general AQCP 
given by J3 = Ji + J2 but Ji 7^ J2 , Eq. (^5)1 implies that 
the dispersion is linear along a direction given by ei = 
a/3( J2 — Ji)i + 3( Ji -I- J2)j and quadratic in a direction 62 
which is perpendicular to ei. Hence i/j^ = 1 and L±^ is the 
length of the system in the ei direction, while v\\ = 1/2 
and i|| is the length of the system in the 62 direction. 

We will now show that the ground state of the model 
in Eq. (jl9p can be written as a product over all k lying in 
half the Brillouin zone. Firstly, the unprimed Majorana 
fermions and must be chosen to have the lower 
eigenvalue EZ of Ht] the corresponding normalized state 



is given by 



1%) 



where e = 



(I/V2) ( a\ 

at + iPt 



^£ bl) 1$) 



"I 



'^k 



(24) 



and 1$) is the vacuum state annihilated by a^, a'p bj: 
and Secondly, the condition Dfi = «&^a^ = 1 for all n 
implies that if we define the Dirac fermion operators c^j = 



6 



(l/2)(a^ — the ground state must be an eigenstate of 
ctcrt with eigenvalue 1 for all n. Hence the state must be 
annihilated by ct for all n; taking the Fourier transform 
of this means that the state must be annihilated by both 
4 = (1/2) (a'i + ib't) and cr = (1/2) (a'- + i6'-) for ah k. 
Hence the normalized state is given by 

|Tj) = (1/V2) ( a'-f + * b'h 1$) (25) 



which ensures that all the points in the rhombus are cov- 
ered uniformly. Once again, we need to avoid the corners 
of the Brillouin zone (i.e, the values and 1 for vi and 
V2), otherwise the fidelity will turn out to be zero. We 
will let vi and V2 go from 1/(2L) to 1 — 1/(2L) in steps 
of 1/L, where L is a large integer. Finally, we must take 
vi + V2 > 1 so as to restrict the integral to half the Bril- 
louin zone. 



for each k. The complete ground state is therefore given 
by the product 



1*) 



n 



m- (26) 



Using Eq. (|26p , we can write the ground state fidelity in 
the form^ 




where 



cos 



and 



with the ± in the superscripts denoting the corresponding 
values with J3 ± 5. One finds 



■K — TV I L 

/L 



■K — -K I L 



-I 



TV/L 



k 



■PI 



Analyzing for small 5 close to the AQCP, we find 
9(52^2 kldk^dk,, 

Jn/L ■ 



InF ; 



27r2 



where R± — 9ky 



/L 
2 



(29) 



(30) 



[^|fc2 — A ± (5)^, and we have only in- 
cluded contributions coming from the low energy modes 
close to the critical modes and extended the limit of in- 
tegrations to 00. 

In subsequent sections, we will investigate the fidelity 
between the two ground states of the model with inter- 
action terms J3 = J^^c — X + S and J3 = J^^c — A — (5, 
respectively, with Ji = J2, i.e., along the vertical line in 
Fig. [5] here A and 5 determine the location in the phase 
diagram. 

Wc will use the simplified equation (PH)) to derive the 
scaling of fidelity analytically. On the other hand, for 
the purpose of numerical analysis of Eq. ((29)) . we will 
parametrize the momenta k^ and ky in terms of two in- 
dependent variables vi and W2 for < wi, W2 < 1, given 
by 



- V2 



1) and ky 



27r 



ivi^V2), (31) 



B. General scaling of fidelity near an AQCP 

We will now proceed to derive a scaling form for the 
fidelity in the thermodynamic limit near a d-dimensional 
generic AQCP in the same spirit as in Ref. [s^ The cor- 
responding scalings in the limit of small system size is 
given in Ref. Is^ . We consider a situation in which the 
correlation length exponent and system size are given by 
1/ = and L = respectively, along m spatial dimen- 
sions, and 1/ = i/j^ and L ~ L±, respectively, along the 
remaining d — m dimensions. We encounter such a case 
with d = 2, m = 1, i^ii = 1/2 and v± = 1 in the two- 
dimensional Kitaev model (point (A) in the phase dia- 
gram) and also near a semi-Dirac band crossing point^. 
We consider the scaling parameter— 

S{X + 6,X-S) 

\n\{M^ + d)\M>^-5))\ 



— lim 

— lim 



N 

\nF{X-S,X + S) 
N ' 



(32) 



where N = L'^L^"' is the system size, A is the distance 
from the AQCP, and A, S are assumed to be positive. We 
propose the scaling ansatz 



S{X + 5,X-S) = 



L-"'L~'"~"'' /((A + (5)L||' II , (A + (S)iy ^, 



{X-6)L](-'\iX~S)LY''^), 



(33) 



where / is a scaling function that is symmetric with re- 
spect to the operation S — )■ —S. Rescaling L||(Lj_) to 

611 {b±) and choosing bn, b± such that (A -I- 6)b 



(A + 6)b]^''^ ^ 1, we get 



i/'^ii 
II 



S{X + S,X-d) 
= (A + (5)''ii™+''^(rf-™)/ [^1, . (34) 

Taking the limit 5/X 0, and expanding f (^^1 T+^/x) ^ 
g{S/X) around (5/A = 0, we arrive at the scaling form 



\nF{X + 6,X-S) 



(35) 
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where we have taken g{0) = g'{x)\x=o- Now let us focus 
on the case A = 0, i.e., we are studymg the fidehty be- 
tween two states at S and —S, respectively, on either side 
of the AQCP. In the limit A = 0, Eq. dM]) shows that 



\nF{S,-S) 



(36) 



Wc note that the above scaling forms arc valid only as 
long as the corresponding exponent of A (see Eq. (|35|) ) 
or 6 (see Eq. ([5^ ) does not exceed 2. Otherwise the 
low-energy singularities associated with the critical point 
become subleading to the quadratic scaling form of per- 
turbation theory, and | lnF| starts varying as A^ (or as 
if A = 0) instead, irrespective of the critical exponents^^. 
Both Eqs. (|55|) and reduce to the scaling presented 



for I'l 



II 



in Ref. [3 

Now wc will 
thermodynamic limit 



consider 



1 the non- 
and choose 



Eq. 

{S « Ll''^^) 

jj_ ^ ^11 ■ ^^^^ limit, a cross-over from a de- 
pendence on (5 to a dependence on the L± takes place in 
the scaling in Eq. p6p which then takes the formal 



\iiF{S,-d) 



—rn 



(37) 



where the (5^ in Eq. ([37]) arises due to perturbation theory. 

In contrary, when S <^ L|| and L|| ^^'^^^ > L~^^"^ ^ we 
get 



\nF{5,-5) 



-5'L 



2j{d-m) - 



- {d—m) 



± 



II 



(38) 



We will now verify the above scaling for the AQCP 
(A) shown in Fig. [5] and determine the fidelity between 
the two ground states at J3 = J^^c + S and J3 = Js^c — S 
with Ji = J2; the system lies in the gapless phase for 
J3 — J^^c ~ S and in the gapped phase for J3 = J^ c + ^■ 
For all numerical studies presented hereafter we have set 
i|| = L± ^ L. 

We use Eq. (|30|) to arrive at the scaling relations 
followed by the quantum fidelity; rescaling kx ^ = 
kxl and ky ^ ky = ky/6, we get 



InF 



27r2 

9^3/2^2 
_^3/2^2 



tt/lVS Jtt/LS 



ky dk^dky 



ky dk^dky 



(39) 



in the limit 6 ^ L ^/"^ = L ^. as expected from Eq. 
([55)1 (see Fig. [3] for numerical verification) . In the above 

Eq. dSll) we have taken i?± = 9k'y^ + (^^k^''^ ± 1 

In the non-thermodynamic limit of J ^ L~^, on the 
other hand, we can use the transformation k^ 



to arrive at the scaling 



lnF« - 



9S^ 
27r2 



fr2 



k^yP+P^ 



(40) 



where P± = 9+{3q^/4±S/ky) . Now, our scaling trans- 
formation suggests q ^ 1/ s/Tj <^ 1 and also the limit 
6 ^ implies 5/ky ^ 1. The above analysis shows 
that for small values of q and ky, which give the domi- 
nant contributions to the integral in Eq. (j40|) . P± are of 
the order of unity. Therefore we get 



InF 



dq 



7t/L ky^^ Jw/Ly^ P+P- 



r°° dk 



3 /2 

Tt/L ky 

^ ~6^L^Xf{J3 = J3,c) - -5^L'>'\ (41) 

which is in complete agreement with our prediction in 
Eq. ((57)) . as shown in Fig.m We note that earlier studies 
of the fidelity susceptibility in the thermodynamic limit 
in the two-dimensional Kitaev model have pointed to the 
same scaling form as in Eq. (|4Tpii2^. 

We reiterate that the study of the scaling of fidelity 
in the thermodynamic limit is closely related to that of 
fidelity per site^^"— . The quantum phase transition at 
J3 = J^ c is associated with a singularity in the double 
derivative of the scaling function S = — lim7v_i.oo luF/N 
given by 



InS 
-djf 



Cln I J3 — J^^d + constant. 



(42) 



where C is a negative constant and hence one observes a 
dip close to the QCP, as shown in Ref. [25l 




FIG. 3: Variation of |lnF[ with 5 as obtained numerically 
at the AQCP in the thermodynamic limit for Ji = J2 = 1, 
L = 1001 and A = 0. \lnF\ varies as 5^/^. Inset: Variation 
of I InFj with L for 5 ^ 0.001 and A = 0. | InFj varies as L^. 
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FIG. 4: Variation of | lnF| with L as obtained numerically at 
the AQCP in the non-thermodynamic limit for Ji = J2 = 1, 
& = 0.00001 and A = 0. | lnF| varies as L^/^ Inset: Variation 
of I ln_F| with 5 as obtained numerically at the AQCP in the 
non-thermodynamic limit for Ji — J2 — 1, L — 1001 and 
A = . I ln_F| varies as 5^. 




0.05 



FIG. 5: Variation of X^^^ I In -F| with A as obtained numerically 
inside the gapless region in the limit (5L S> 1, for Ji = J2 = 1, 
S = 0.001 and L = 3001. |lnF| varies as A^^/^lnA. Inset: 
(a) Variation of | lnF| with L as obtained numerically inside 
the gapless region for Ji = J2 = 1, S = 0.001 and A = 0.01. 
I ln_F| varies as L^. (b) Variation of | ln_F| with 5 as obtained 
numerically inside the gapless region for Ji = J2 — 1, L = 
5001 and A = 0.01. I InFl varies as S^. 



C. Fidelity inside the gapless region: 

In this section, wc consider the situation when both the 
states under consideration lie inside the gapless region of 
the phase diagram (Fig. [2]) along the dashed vertical line 
with A = J3 - J3,c » <5 > and A > L-'^^"^ = L^^. 
To calculate quantum fidelity, we numerically integrate 
Eq. ([29)) and arrive at the scaling relation 



Ini^ - -^2^2^^-1/2 i^-^x 



(43) 



in the limit SL^^'^-^ = (5L ^ 1. In Figs. we present 
the numerical results which clearly support the above 



1. 1x10 

lln Fl/L- 




7x10 



1000 



2000 



3000 



FIG. 6: Variation of | ln_F|/L^ with L as obtained numerically 
inside the gapless region in the limit 5L <^ 1, for Ji = J2 — 1, 
6 = 0.00001 and A = 0.005. I Inf I varies as InL. 



scaling prediction. 

Close to the AQCP, one can provide an analytical veri- 
fication of P5|) using Eq. ([5(1)) with R± as defined before. 
Using the transformations — kx/VX., ky = ky/X^ Eq. 
(PH]) can be rewritten as 



00 poo 



hiF 

'tt/LA Jn/LV\ 

6'L'\-'^'giL,X), 



ky dk^dky 



9ky' + ( 



(44) 



in the limit X ^ S when S appearing in the integrand 
can be ignored. The function g{L, A) is found to scale 
as g{L, A) ^ In A by numerical investigations of Eq. (|29p . 
We interpret this logarithmic behavior in (|43|) as a signa- 
ture of the system being in the gapless region. The scal- 
ing |lnF| - A-i/2lnA can be understood noting that 
A denotes the distance from the AQCP; the power-law 
scaling A^^/2 foUows from the generic scaling in Eq. (|35p . 
while the gapless nature of the phase diagram is encoded 
in the additional logarithmic correction. 

On the other hand, in the limit 6L ^ 1, again with 
A ^ L^^, a similar analysis of Eq. ([29| leads to the 
scaling 



InF 5^L^X~^^^\nL\nX, 



(45) 



as shown in Fig. \6\ we therefore find an additional In L 
correction in comparison to the scaling in (|43p . Interest- 
ingly, it can be shown that there exists another cross-over 
at A ^ L~^^'^^ = L~^, when the system size dependence 
changes to InF ^ U'l'^. This is expected as the system 
approaches the vicinity of AQCP where the scaling Eq. 
((37)) is apphcable. 

A few comments are necessary at this point. Our anal- 
ysis points to a cross-over from \vlF ^ for JL^/"^ = 
8L~>\lo\^F ^ L'^XnL for bL < 1 with A > L"!/'^^ 
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in both the cases (see Eqs. (43) and (45), above). This 
apparently suggests that even in the gapless phase, we 
see a cross-over around 5L ^ 1, which resembles a ther- 
modynamic to non-thermodynamic cross-over in fidelity, 
as observed in Ref. [3^ though scaling with S remains 
the same in the present case. It also appears that the 
crossover occurs as 8L^^^^ ~ 1,, which suggests that the 
AQCP may play the role of a dominant critical point in 
its vicinity in the gapless phase. 



D. Fidelity inside the gapped phase: 

For A < 0, both the states are in the gapped phase 
along the vertical line of Fig. [2l and choosing A ^ L^^ 
and A 3> (5 (i.e., non-thermodynamic limit), one finds 
numerically 

which matches exactly with non-thermodynamic result 
for fidelity susceptibility^^. As discussed above, we en- 
counter a cross-over to | lni^| ^ L^^"^ for A ^ L~^. 




1400 2000 L 3000 

FIG. 8: Variation of | ln_F| with L in the non-thermodynamic 
limit, as obtained numerically for Js = Ji + J2, Ji = SJ2 = 3, 
A = and 5 = 0.00001. | ln_F| varies as L^^^ in this regime. 
Inset: Variation of |lnF| with S in the non-thermodynamic 
limit, as obtained numerically for J3 = Ji + J2, Ji = 3J2, 
A = and L = 1001. | \nF\ scales quadratically with S. 



IV. CALCULATING FIDELITY IN THE 
KITAEV MODEL USING ROTATION OF SPINS 



E. Observations for Ji 7^ J2 

Our attention so far has been concentrated on the case 
Ji = J2- However, all the points on the critical line J3 = 
Ji + J2 correspond to an AQCP, regardless of whether 
Ji = J2 or not. In Figs. [7] - [H we have presented our 
numerical results with Ji = 3 J2 which clearly shows the 
effect of the AQCP. 




0.002 5 0.01 



FIG. 7: Variation of | In F\ with 5 in the thermodynamic hmit, 
as obtained numerically for J3 = Ji + J2, Ji = 3J2 = 3, 
L = 3001 and A = 0. | ln_F| scales as S^^^ . Inset: Variation of 
I In _F| with L in the thermodynamic limit, as obtained numer- 
ically for J3 = Ji -f J2, Ji = 3J2 = 3, A = and S = 0.001. 
I ln_F| shows a quadratic scaling with L in this limit. 



In this section, we will compute the overlap between 
two ground state wave functions of the Kitaev model with 
each spin rotated about some axis by an angle 77 and 77 4- 
drj, respectively; we note that a similar method has been 
used to calculate the geometric phase close to a QC P^°'^^ . 
Let us recall the complete ground state given by the prod- 
uct form in Eq. (^5]) . To compute the fidelity, let us 
introduce a family of Hamiltonians generated by rotat- 
ing each spin by an angle 77 about the z directioniSi^i^, 
i.e., H{r]) = g^Hg^ with .g,, = Hn exp(i?7fT£/2); this 
unitary transformation leaves the energy spectrum in 
Eq. unaltered. Under this rotation, the Pauli spin 
matrices transform as (t5 — > cos 77 crS — sin 77 and 
crti — >• cos 77 crH + sin 77 crS. Hence the Majorana fermions 
transform to 

0^5(77) = cos 77 an + sin 77 a^, 
O'ni'n) = cos 77 - sin 77 a^;, 
bfiiv) = cos 7763 - sin 77 6^, 
6^-5(77) = cos 77 b\i + sin 77 hn, (47) 

with similar expressions for aj:{ij), '^\{v)- etc. The ground 
state of H{ri), denoted by 1^* (»?)), is therefore given by 
an expression similar to Eq. (P^ . with at, a'^, 6t, b'~ 

being replaced by at (77), a'^iv): ^^('')- 

We now find that the overlap between the ground 
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states for two different values of ry is given by 



where sin6'r = 



"I + ^1 



(48) 



(Note that the overlap is unity for both 772 = fji f^^d 
= ?7i + tt) . Eq. (|3S|) implies that 



ln(*(0)|«'(dr/)) 



1 - sin6lg) 



^'^''^^ i d^fc (l-sin0,,), 



(49) 



up to order (dr?)^, where A = 47r^/(3A/3) denotes the 
area of half the Brillouin zone over which the integration 
is carried out in the second equation in (We recall 

that the number of k points in half the Brillouin zone is 
given by i/4). 

The above expression for the fidelity shows that there 
is no term of first order in drj in the present case; hence 
the geometric phase is zero. The coefficient of the second 
order term, dr]"^, yields the fidelity susceptibility. Note 
that this is proportional to L and does not exhibit any 
non-analytic behavior as a function of the couplings Ji, 
J2 and J3. This is not surprising; a rotation of all the 
spins is simply given by a unitary transformation, and 
the system does not cross a QCP as a result of such a 
transformation. 



one-dimensional system of massive Dirac fermions with 
and without interactions and in the Kitaev model on the 
two-dimensional honeycomb lattice. The behavior of the 
fidelity in the one-dimensional Dirac system agrees with 
the general scaling predictions made earlier—. We have 
also derived general scaling relations for the fidelity close 
to an AQCP and have verified our predictions by using 
the AQCP present in the Kitaev model. Moreover, we 
observe an additional logarithmic correction (in the lin- 
ear dimension L of the system) in the scaling form of the 
fidelity inside the gapless phase of the two-dimensional 
Kitaev model when SL^^'^^ <^ 1. Our numerical stud- 
ies apparently indicates a crossover in scaling around 
SL^^"^ - 1. Finally we have considered a rotation of all 
the spins in the Kitaev model by an angle 77 about z-axis 
and calculated the fidelity between two ground states cor- 
responding to two different values of rj. Wc have shown 
that the geometric phase is absent and the fidelity does 
not show any singularity because no QCP is crossed when 
such rotations are performed. 

Acknowledgments 



AD and VM acknowledge Ayoti Patra for collaboration 
in related works. AD acknowledges CSIR, New Delhi, for 
financial support and DS acknowledges DST, India for 
Project SR/S2/JCB-44/2010. 



V. CONCLUSION 

We have studied the ground state fidelity in both the 
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